Hydrogen Dissociation and Diffusion on Ni and Ti -doped Mg(0001) Surfaces 
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It is well known, both theoretically and experimentally, that alloying MgH2 with transition ele- 
ments can significantly improve the thermodynamic and kinetic properties for H2 desorption, as well 
as the H2 intake by Mg bulk. Here we present a density functional theory investigation of hydrogen 
dissociation and surface diffusion over Ni-doped surface, and compare the findings to previously 
investigated Ti-doped Mg(0001) and pure Mg(0001) surfaces. Our results show that the energy 
barrier for hydrogen dissociation on the pure Mg(0001) surface is high, while it is small/null when 
Ni/Ti are added to the surface as dopants. We find that the binding energy of the two H atoms near 
the dissociation site is high on Ti, effectively impeding diffusion away from the Ti site. By contrast, 
we find that on Ni the energy barrier for diffusion is much reduced. Therefore, although both Ti and 
Ni promote H2 dissociation, only Ni appears to be a good catalyst for Mg hydrogenation, allowing 
diffusion away from the catalytic sites. Experimental results corroborate these theoretical findings, 
i.e. faster hydrogenation of the Ni doped Mg sample as opposed to the reference Mg or Ti doped 
Mg. 



PACS numbers: 

I. INTRODUCTION 



Safe and efficient hydrogen storage is one of the biggest 
barriers to the more wide spread usage of hydrogen as an 
energy carrier or fuel. Currently, commercial solutions 
are based on liquid or compressed gas storage methods, 
which are inefficient and have safety issues. Alterna- 
tive storage methods include metal hydrides, which are 
formed by the interaction between a suitable metal and 
hydrogen. The relatively strong metal-hydrogen bonds 
provide an intrinsically safe storage medium. The re- 
lease of hydrogen from the hydride is then achieved by 
heating the material above a certain decomposition tem- 
perature. There are a large number of metals in nature 
that form hydrides, however, only the lighter ones are 
thought to be suitable candidates for mobile hydrogen 
storage purposes (see [5] for an overview). Beside being 
light weight , a hydride will need to have good cyclability 
(several hundred times with little loss of performance), 
fast adsorption/desorption kinetics (the hydride should 
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form/decompose on a time scale of minutes) and low 
decomposition temperature (ideally between 20 and 100 
Celsius). 

Magnesium is a good case study due to its lightweight, 
low cost, cyclability and the high H storage capacity of 
7.6% by weight once the hydride MgPi2 is formed [3]. 
However, its commercial application is still on hold for 
practical issues due to low H absorption/desorption ki- 
netics and high working temperatures [I]. The strong 
bond between Mg and hydrogen provides Mgfi2 with 
high thermodynamic stability, which has an enthalpy of 
formation of about -76 kJ/mol jS], and a decomposition 
temperature of more than 300 Celsius [6 . The slow ki- 
netics may be explained by the high energy barrier which 
needs to be overcome (see for example [7]) to dissociate 
the H2 molecule due to the tendency of Mg to repel the s- 
electron of H because of the Pauli exclusion principle [8 . 

A step forward in improving hydrogen reaction kinetics 
has been achieved by the mechanical ball milling of Mgfi2 
with transition elements (see [5] and references therein). 
The hydrogen storage properties of mechanically milled 
powders improves because of the reduced powder size (see 
for example [TTjl ITTl HB] and references therein) , which 
shorten the diffusion distance of H into bulk Mg for the 
formation of the hydride. There are many experimental 
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and theoretical papers in the literature showing that the 
hydriding properties of MgH2 are further enhanced by 
the addition of traces of transition metals which act as 
a catalyst (see for example [T3J HH [TSJ [IDl [T7] and ref- 
erences therein). In particular, alloying Mg with Ni can 
slightly improve the thermodynamic properties of the hy- 
dride by favouring H2 adsorption/dissociation and con- 
sequent atomic hydrogen absorption/desorption due to a 
weakening of the bonding between Mg and H atoms (see 
for example [TT1 IT31 HF1 ITfl ITB] and references therein). 

While from theoretical calculations the destabilization 
effect of Ni on MgH 2 appears second only to Cu (see [IB]), 
experimentally Ni shows the highest kinetics, with Cu 
falling behind. As suggested by Shang et al. [13], Cu re- 
sults are disadvantageous for H desorption probably be- 
cause of the formation of a MgCu2 compound. Recently, 
a new method of chemical fluid deposition in supercrit- 
ical fluids has been used on metal hydrides [19]. Even 
sparser literature exists for the activation barrier of hy- 
drogen dissociation on a transition metal doped surface, 
which includes only the theoretical calculations made by 
Du et al. [TJ EQl E] within DFT (RPBE) for both the 
pure Mg(0001), and Ti and Pd incorporated Mg surfaces. 
Their results show that the dissociation barrier of hydro- 
gen on the Ti doped Mg surface is greatly reduced (in 
fact, there is no barrier at all) due to the strong interac- 
tion between the hydrogen s orbital and the Ti d orbital, 
however, strong binding of the two H atoms near the Ti 
site prevents easy diffusion, reducing therefore the effi- 
cacy of the catalyst for Mg hydrogenation 1 . Palladium 
doping appears to both lower the dissociation barrier and 
the diffusion barrier, suggesting a better catalytic activ- 
ity. Their findings are consistent with the experimentally 
observed trend of generally improved hydrogen absorp- 
tion kinetics when Mg surfaces are doped with transition 
metals, as previously mentioned. 

To our knowledge, so far there is no published theo- 
retical investigation of H2 dissociation and correspond- 
ing activation barrier on a Ni-incorporated Mg surface, 
nor a systematic investigation of the catalytic effect of 
other transition metal dopants apart from the above men- 
tioned studies on the Ti-doped and Pd-doped Mg sur- 
faces presented by Du et al. [TJ [2DJ [21]. There are a 
few theoretical papers about the dissociation of molec- 
ular hydrogen on a pure Mg surface where the corre- 
sponding activation barrier has been effectively calcu- 
lated. These investigations were based on a jellium model 
and potential energy surface (PES) calculations within 
density functional theory (DFT) with the local density 
approximation (LDA) or generalised gradient corrections 

(rpbe) [2U EH ESI • 

For the purpose of a larger scale investigation, we have 
performed DFT calculations for hydrogen dissociation 
and diffusion on a Ni-doped Mg surface, accompanied 
by analogous calculations on a Ti-doped Mg surface for 
a consistency check with the recently reported theoreti- 
cal values. This study should be regarded as a first step 
in order to build up a global picture of the dissociative 



chemisorption of hydrogen when doping the Mg surface 
with different transition metals. The main purpose of 
this article is to try to understand the observed large 
enhancement of the kinetics of hydrogen adsorption by 
Mg when it is doped with a small quantity of Ni, but 
not when it is doped with Ti. The computational results 
are supported by experimental data where a Ni doped 
Mg sample is hydrogenated substantially faster than the 
reference Mg or Ti doped Mg. 

II. COMPUTATIONAL METHOD 

DFT calculations were performed with the ab-initio 
simulation package VASP [27] using the projector aug- 
mented wave (PAW) method (28J [29] and the PBE 
exchange-correlation functional [3D]. An efficient charge 
density extrapolation was used to speed up the calcula- 
tions |31j . A plane- wave basis set was used to expand the 
electronic wave-functions with a plane-wave energy cut- 
off of 270 eV, which guarantees convergence of adsorption 
energies within 1 meV. For completeness, Mg bulk pa- 
rameters were also calculated using the LDA functional. 
Monkhorst-Pack k-points were used to sample the Bril- 
louin zone [32 . A smearing function of the Methfessel- 
Paxton (MP) type (product of a Gaussian times a nth- 
order Hermite polynomial) |33j was used throughout. 
Figs. [T][3]-|9]andp]have been made using the XCRYS- 
DEN software [34 ■ The exact values of the various pa- 
rameters used in the calculations will be reported below 
in the relevant sections. 

Activation energies have been calculated using the 
nudged elastic band (NEB) method [35]. A sufficient 
number of replicas has to be used in order to predict ac- 
curately a minimum energy path (MEP), for most cases 
we repeated the calculations with a different number of 
replicas until convergence of the activation energy and 
main features of the MEP were observed. The total num- 
ber of images actually used in each case is reported where 
relevant in the following sections. 

III. THEORETICAL RESULTS 

A. Bulk Mg, Ti and Ni, and the Mg(OOOl) surface 

Magnesium bulk crystal at ambient conditions has the 
hexagonal closed packed (hep) structure. Several prelim- 
inary tests were first carried out using the PBE version 
of the PAW potential of magnesium, which only includes 
the 3s 2 electrons in valence and has a core radius of 1.1 A. 
These included: the energy dependence on the c/a ratio 
for different k-points meshes, from a minimum of 56 to 
a maximum of 880 k-points in the irreducible wedge of 
the Brillouin zone (IBZ); different values of n for the MP 
smearing functions and different smearing widths; and 
different plane-wave cutoffs. To calculate the bulk struc- 
tural properties of Mg, energy versus volume curves were 
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fitted to a Birch- Murnaghan equation of state [35] . We 
found that with a 18x18x12 k-point mesh (259 points in 
the IBZ), n = 1 and a smearing width of 0.2 eV, the 
zero pressure equilibrium volume Vq and bulk modulus 
Bq of bulk Mg were converged to within 0.2% and 0.3% 
respectively. Similar convergence results were obtained 
with the standard LDA potentials, which also has only 
the 3s 2 electrons in valence and a core radius of 1.1 A. 
Results are summarised in Table |T] together with results 
from previous theoretical calculations. The well known 
trend of LDA to overestimate the bulk modulus and un- 
derestimate the lattice parameter |37j is apparent. 

Finally, we have tested PBE and LDA PAW potentials 
with 2p 6 and 3s 2 electrons in valence. These potentials 
still have core radii of 1.1 A but require an higher en- 
ergy cut-off value of 350 eV. The structural parameters 
obtained with these potentials are essentially identical to 
those obtained with the previous potentials. Therefore, 
in the rest of the work we only used the standard Mg 
potentials. 

From our PBE calculations (see Table [I]), we derive 
a lattice constant of 3.19 A, in error of just 0.6% with 
respect to the experimental value [39 . The zero pressure 
bulk modulus Bq is 36.8 GPa, and the value of c/a at the 
equilibrium volume is 1.621, both in very good agreement 
with the experimental values (we note however that these 
calculations do not include room temperature thermal 
expansion, which are present in the experimental data). 

Titanium is also hep crystal. We used the standard 
version of the PBE PAW pseudopotential for Ti, which 
has a core radius of 1.2 A, and we used a 18x18x12 k- 
point grid. The resulting values for the structural pa- 
rameters were clq — 2.923 A and Bq = 120 GPa, and the 
value of c/a at the equilibrium volume was 1.583, in good 
agreement with those previously found with theoretical 
and experimental investigations (see Refs. I4T1 B51 I4"9l and 
references therein) . 

To study bulk Ni we also used the standard version 
of the PBE PAW potential, which has a core radius of 
1.1 A. Ni bulk has a face centred cubic crystal structure, 
with a small magnetic moment of 0.61 hb under ambi- 
ent conditions [42 , so we performed spin polarised to- 
tal energy calculations. The calculations were performed 
with a 13x13x13 grid of k-points. We found a lattice 
parameter a = 3.524 A and a bulk modulus Bq = 194 
GPa, which compare well with the experimental data of 
3.524 A and 186 GPa respectively [42J, [50] . The zero pres- 
sure magnetic moment is 0.63 [is, which is also close to 
the experimental value of 0.61 hb- The values we found 
for do, Bq and /ib are in agreement with those from other 
GGA and PBE calculations [EEJ [52j [53] . 

Surfaces have been modelled using periodic slabs, with 
several atomic layers (from 3 to 13) and a large vacuum 
thickness (5-18 A), defined as the distance between two 
opposite facing surfaces. We used a 18x18x1 k-point grid 
for the lxl surface primitive cell. The positions of the 
atoms in the three topmost layers were allowed to relax, 
while the rest were kept at the bulk interatomic distances. 



Good convergence in the calculated surface energies and 
relaxations of the topmost atomic layers was achieved 
with five layer slabs (corresponding to a slab thickness 
of about 13 A) and a vacuum region thickness of about 
10 A. We found that the topmost layer has an inward re- 
laxation of about 1.4 %, in good agreement with the in- 
ferred experimental zero temperature value of 1.7% [54] . 

We found that with 5 atomic layers the surface energy 
is converged to within 2 meV to the value 0.30 eV/atom. 
This compares well with the experimental findings which 
are in the range 0.28-0.33 eV/atom [4"5ll4*6] , 



B. H2 dissociation and H diffusion on the pure Mg 
and the Ti-doped Mg surfaces 

Hydrogen adsorption energies on the Mg(0001) surface 
were determined at low coverage in four possible sites: 
top, bridge, hollow-fcc and hollow-hep (see Fig. [I]). These 
adsorption energies are defined as E ads (H)— i? slab (MgH)- 
[£ s i ab (Mg) + 1/2 E{R 2 )], where £ s i ab (MgH) is the en- 
ergy of the slab with one H adsorbed on the surface, 
-Esiab(Mg) is the energy of the bare slab and E(H2) the 
energy of the isolated hydrogen molecule, calculated by 
placing the H2 molecule in a large cubic box of sides 
13.5 A. 

Calculations have been performed on 2x2 (correspond- 
ing to 0.25 ML coverage, ML=monolayer) and checked 
against results obtained from 3x3 (corresponding to 0.11 
ML coverage) surface unit cells, with differences between 
the two sets of calculations of less than 0.01 eV, thus 
implying that the results effectively correspond to those 
for an isolated H2 molecule. The two sets of calculations 
have been performed with equivalent grids of k-points, 
9x9x1 and 6x6x1 for the 2x2 and the 3x3 surface unit 
cells respectively. 

The values for the adsorption energies of atomic hydro- 
gen in different adsorption sites on the pure Mg surface 
are reported in Table |TTJ These compare well with pre- 
vious theoretical results 25j. It is clear that there is a 
strong preference for the hollow sites, with a small pref- 
erence for the fee hollow site. 

We performed NEB calculations for H2 dissociation 
over two possible sites (bridge and top) . These have been 
accompanied by careful tests on supercell size, number of 
layers in the slab and number of replicas in the NEB cal- 
culation to obtain the minimum energy path and the ac- 
tivation barrier. We found that a 2x2 supercell, 5 layers 
and 5 replicas are enough to obtain activation energies 
converged to within 0.02 eV. The first MEP is rather 
featureless (IS — > TS — > FS), and it is well approximated 
also by 5 replicas, although in Fig. [2]we report the MEP 
obtained with 17 replicas. 

Of the two sites investigated, we found that H 2 prefers 
to dissociate over a bridge site (see Fig. [3| with an acti- 
vation energy of about 0.87 eV (about 0.6 eV lower than 
that obtained for dissociation over a top site), in agree- 
ment with previous DFT calculations [2"0"1 [Ml [231 12"6] and 
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TABLE I: Bulk and surface properties of pure Mg. 



ajA) c/a ^ell/atom) ^ k (GPa) dkp 

This work 3.19(3.13)" 1.621(1.621)" 22.85(21.59)" 36.8(40.5)" 3.9(4.3)" 

Other calculations 3.19 6 , 3.18(3.13) c , 1.615(1.616) c , 22.97(21.66) d 35.5(40.2) c , 4.0(4.1) d 

3.20(3.14) d 1. 624(1. 622) d 36.0(40.1) ti 

[Expt.] [3.21] e [1.624] c [23.24] / [35.4] e , [36.8 ± 3.0] 9 [4.3 ± 0.4]'- 

E coh {eV) -1.50", -1.50(-1.78) c 

[Expt.] -1.51' 1 

E surf (ev/atom) 0.30", 0.30(0.35) c , 0.34 s , 0.32 3 

[Expt.] 0.28*, 0.33' 



"Reported values are those from PAW PBE(LDA) calculations 
which do not include room temperature thermal expansion. 
6 Ref. [201 

c Ref. [371 from DFT GGA(LDA) calculations. 
d Ref.[38lfrom PAW GGA(LDA) calculations. 
e Ref. [39] 
•^Ref. \M 
9 Ref. EH 
^Ref. 1321 

'Ref . B3| from ab initio LDA calculations. 
3 Ref . R4l from ab initio LDA calculations. 
fc Ref. 05] 
'Ref. \M 



TABLE II: Hydrogen adsorption energies (E a( ; s ) in different 
adsorption sites on the pure Mg surface, for the 2x2 and the 
3x3 surface unit cells. 



TABLE III: Activation energy (E a ) for hydrogen dissociation 
on the pure Mg, Ni-doped and Ti-doped Mg surfaces. 



Ads. sites 


E ads (eV) 


E ads (eV) 




(2x2) 


(3x3) 


Top 


0.75 


0.74 


bridge 


0.12 


0.13 


hollow (hep) 


-0.03 


-0.03 


hollow (fee) 


-0.05 


-0.04 



experimental findings [7j, as reported in Table III 

The small difference between our findings and those of 
Vegge [T7] and Du et al. [5D] are due to their use of RPBE 
instead of PBE, and different k-point meshes. 

We then performed a second NEB calculation to ob- 
tain the MEP for the diffusion of one of the H atoms 
on the surface from one fee to a second fee site (FS — > 
TS2 -> LS -> TS3 -> FS2; see Fig. H. This MEP (cal- 
culated with 17 replicas) is also displayed in Fig. [2| as 
a continuation of the dissociation MEP, and shows that 
the highest energy barrier for surface diffusion is only 
^0.18 eV, which agrees very well with the calculations 
of Du et al. [1 . This low energy barrier clearly indicate 
fast diffusion even at room temperature. 

Before repeating the calculations on the Ti-doped sur- 
face, we tested all four possible sites for H adsorption 
after dissociation (see Fig. [5]), and we found that atomic 
hydrogen prefers to adsorb into two of the possible three 
hollow-fcc sites around the Ti atom. The dissociation 



E a (pure Mg) 
[Expt.] 

E a (Ni-doped Mg) 
E a (Ti-doped Mg) 



0.87", 0A b - c , 0.5 d ' c , 1.15^ 1.05 9 , 0.95 fe 
1.0' 

0.06" 
null", negligible 9 



"This work. 

'Ref. 22 for a jellium system. 

c Ref.[25] from DFT LDA calculations ans PES. This lower value 
as compared to other calculations is explained as due to the well 
known LDA overbinding. 

d Ref. 23 for a jellium system. 

e Ref . [24l for a jellium system and PES. 

^Ref. [261 from DFT RPBE. 

9Ref. [20j from DFT PAW RPBE calculations (see also discussion 
in main text). 
h Ref. 55 from PES calculations. 
! Ref. H 



activation barrier was calculated using 9 and 17 replicas, 
with 9 being enough to display the main features of the 
MEP (IS -> FS; see Fig. [§, although in Fig. [| we dis- 
play the results obtained with 17 replicas. Our findings 
are very similar to the previous results of Du et al. [T|l2"0]. 
i.e. there is no barrier for hydrogen dissociation on a Ti- 
doped Mg surface, and a barrier of almost 0.8 eV for 
diffusing away from the Ti sites (FS TS2 FS2; see 
Fig.|7|, which therefore becomes the rate limiting step in 
the reaction [T]. 
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H2 dissociation and diffusion on a Ni-doped Mg 
surface 



Having benchmarked our calculations on the pure Mg 
and the Ti-doped Mg surfaces, we now come to the main 
purpose of the paper, which is to study the effect of Ni 
doping of the Mg(0001) surface on the activation barriers 
for H 2 dissociation and diffusion on the surface. 

On the Mg(0001) surface, we found that Ni is non- 
magnetic, so all calculations have been performed with- 
out including spin-polarisation. 

After dissociating on top of a Ni atom, the two H atoms 
can adsorb into four different hollow sites, as shown in 
Fig. [5j The most stable final state is found to be the one 
where the H atoms adsorb into two nearby hollow-hep 
sites (see Fig. [5] bottom-right corner. We also found that 
the configuration on the bottom-left corner was unstable, 
with the hydrogen atoms repelled by the Nickel atom and 
squeezed between nearby Mg atoms) . Figure [8] shows 
the dissociation of the hydrogen molecule over the Ni 
atom as viewed from side (top panel figures) and top 
(bottom panel figures) positions respectively at the IS, 
TS and FS. Note that on the Ni-doped Mg surface the 
molecule at the TS is much higher than on the pure Mg 
(0001) surface (see respectively Figs. [8] and [3]), being at 
~ 2 Aand ~ 1 A over the two surfaces respectively. 

NEB calculations were run with different numbers of 
replicas, and we found that 9 replicas are enough for a 
precise calculation of the energy barrier to within 1 meV. 
The resulting activation barrier for H2 dissociation on a 
Ni-incorporated Mg surface is only 0.06 eV, against 0.87 
eV found for the pure Mg surface. In Fig. [2] we display 
the MEP obtained from a calculation with 17 replicas (IS 
-> TS -> FS). 

The NEB diffusion calculation was also performed with 
17 replicas (FS -> TS2 -> FS2; see Fig. ^ and shows 
an energy barrier of only 0.27 eV, which is only slighlty 
higher than the diffusion barrier on the pure Mg surface, 
and would also allow fast diffusion even at room tem- 
perature. We note that this barrier is similar to the one 
found on the Pd-doped surface by Du at al. [T], although 
in that case the rate limiting step is the dissociation of 
the H2 molecule with an energy barrier of 0.305 eV. 

This suggests that Ni should be an even better catalyst 
than Pd for the hydrogenation of Mg. 

As a final note we would like to point out an inter- 
esting analogy with H2 dissociation on pure transition 
metal surfaces. In particular, on the pure Ni(lll) sur- 
face Kresse [56] calculated an energy barrier of only 0.015 
eV using DFT PAW GGA. This is similar to our value of 
0.06 eV on the Ni-doped Mg surface, however, we note 
that when the same 4x4x1 k-point sampling grid is used 
we find an energy barrier of 0.014 eV on the Ni-doped 
Mg surface, which is therefore very close to the value 
found by Kresse [55]. These calculations are also con- 
sistent with potential energy calculations of Arboleda et 
al. [55], also performed with a 4x4x1 k-point grid. The 
small barrier for hydrogen dissociation on Ni(lll) is also 



confirmed by experiments |57j). 

Analogously, the behaviour of H2 dissociation over the 
Ti-doped Mg surface appears to be similar to that ob- 
tained on the pure Ti(lll) surface: the null activation 
barrier we find with a smaller 4x4x1 grid compares with 
theoretical results found over a Ti (0001) surface with an 
analogous grid [55, 58 . In other words, this seems to sug- 
gest that the value of the activation barrier for hydrogen 
dissociation over a transition-metal doped Mg surface is 
similar to the activation barrier for H2 dissociation over 
the corresponding pure transition-metal surface. 



D. Electronic structure 

To study the electronic properties of the system, we 
projected the electronic density of states onto spherical 
harmonics functions of type s, p and d, centred on Mg, 
Ni, Ti and H atoms. It is well known that the catalytic 
reactivity of a surface is correlated to the position of the 
ri-band (i.e., in this case the projection of the electronic 
density of states onto ri type spherical harmonics) with 
respect to the Fermi energy Ef. In particular, it was 
shown by Hammer and Norskov [5] (see also [59]) that a 
convenient parameter to monitor the catalytic reactivity 
is the first energy moment of the ri-band, or ri-band cen- 
tre, defined as E d = J E ^ dE(E - E f )p d (E), where p d (E) 
is the density of states projected onto spherical harmonic 
of type d centred on some specified atom, and E is some 
cutoff energy which we chose to be at 7 eV above the 
Fermi energy. Then, if the centre of the band is close to 
Ef it follows that there are many d electrons available 
for donation, as well as a significant number of empty ri- 
levels available for back-donation, and the results of this 
is that the system is very reactive. The d electrons of Ni 
on Mg(0001) form a band which is relatively close to the 
Fermi energy and for this reason the system is very reac- 
tive. By contrast, Mg has no d electrons (although in the 
solid state a projection onto d type spherical harmonics 
is not zero) , and therefore its reactivity is much reduced 
by comparison. 

Ni is a late transition metal with almost all the d or- 
bitals filled with electrons, by contrast Ti only has 2 elec- 
trons in the d orbitals. It is therefore clear that the posi- 
tion of the ri-band centre will be much higher in Ti than 
in Ni, which explains the higher reactivity of Ti. 

Our calculated values for E d on the Ni/Ti-doped Mg 
surfaces are -0.79 eV and +1.08 eV for Ni and Ti respec- 



tively (see Table IV I 



Using DFT RPBE, Vegge et al. [T7] calculated the ri- 
band center position for MgTM (TM = transition metal) 
alloys, allowing an expansion of the alloy lattice to acco- 
modate the hydrogen atoms. They found -0.82 and +0.48 
for TM=Ni,Ti respectively. Although the value we find 
for Ti is much larger, the same trend is observed in the 
case of our Ni/Ti-doped Mg surfaces. 

Figures |10]|T2"] show the projected density of states 
(PDOS) found for the pure Mg, Ni-doped and Ti-doped 
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TABLE IV: The d-band center position with respect to the Fermi energy (Ed — Ej), H a peak shift between the initial 
and transition state (H S TS_/S ), activation barrier (E a ) and energy difference between the final and initial state (E FS ~ /S ) for 
hydrogen dissociation on the pure Mg surface as opposed to the Ni/Ti-doped Mg surfaces. 



Surface 


E db - E F (eV) 


H S TS "' S (eV) 


E n (eV) 


E FS-IS (eV) 


Mg pure 




-1.43 


0.87 


-0.04 


Ni-doped Mg 


-0.79 


-0.77 


0.06 


-0.66 


Ti-doped Mg 


+1.08 




null 


-1.34 



12 



10 



Mg surfaces respectively. The PDOS are given for a num- 
ber of configurations along the MEP: the initial state 
(IS), the transition state (TS), the replica just after the 
transition state (TS+1) and the final state (FS). For the 
Ti-doped Mg surface, the PDOS are given for IS and FS 
only since there is no transition state in the dissociation 
process. For somplicity of notation, we call here the tran- 
sition state and the final state simply TS and FS, as we 
only refer to the part of the MEP which deals with the 
dissociation of the H2 molecule. 

In the IS the hydrogen molecule is still far from the 
metal surface and there is no overlap between the H 2 
molecular orbitals and the orbitals of the metal surface. 
At the transition state, instead, when gaseous hydrogen 
has started dissociating over the surface, there is clear 
interaction between the H s orbital and the Mg s and p 
orbitals on the pure Mg surface (see Fig. 10 top-right 
corner). On the Ni-doped surface the overlap appears 
to be non-zero with all the Ni s, p and d orbitals (see 
Fig. [Til top-right corner). In the final state it is evident 
that the magnitude of the Mg p electron peaks below the 
Fermi level are increased in the Ni/Ti-doped surfaces (see 
respectively Fig. |11[ bottom-right corner, and Fig. 
right) with respect to the pure Mg surface (see Fig 
bottom-right comer). 

Interestingly, we note that there appears to be a clear 
negative shift of the position of the hydrogen s orbital in 
going from the initial state to the transition state, which 
is more pronounced for the pure Mg surface as opposed 
to the Ni-doped Mg surface. 

Besides the <i-band center positions, in Table |IV| we 
also report the corresponding activation barriers (E a ), 
the energy difference between the initial and final states 
(E FS ~~ IS ) for hydrogen dissociation, and the H s peak 
shift between the initial and transition states (H- S TS ~ IS ). 
The correlation between the position of the <i-band and 
the height of the activation barrier is evident, as well 
as the correlation with E FS ~ IS , i.e., the <i-band center 
position is smaller for larger values of the former and 
smaller values of the latter. 

Furthermore, from the results obtained for the pure 
Mg and Ni-doped Mg surfaces, another interesting corre- 
lation emerges. In fact, as shown in Table [rv] it appears 
that R S TS ~ IS correlates with both E a and E FS ~ IS , i.e., 
it is smaller for smaller values of the former and for 
larger values of the latter, following a reversed trend with 
respect to that noticed for the d-band center position. 
In other words, this means that the shift of the hydro- 



gen s orbital between initial state and transition state is 
larger when the bond between the dopant and H atoms 
is weaker. 



E. Charge distribution 

To conclude our analysis we decided to have a look at 
the charge distributions in the systems as the dissociation 
processes take place on the pure Mg and metal-doped Mg 
surfaces. 

To do this, we calculated the total charge at each 
step of the MEP, and for convenience, we subtracted 
the charge densities obtained from calculations which in- 
cluded only the substrate and only the H2 fragments re- 
spectively, with the atoms in exactly the same positions. 
This charge difference obviously integrates to zero, and 
has the advantage of showing point by point where the 
charge is being transferred to. The analysis reveals some 
interesting effects. In particular, on the pure Mg(0001) 
surface we find that at the transition state there is a sig- 
nificant charge transfer from the Mg substrate to the H 
atoms (see Fig. 13 - left). This extra charge fills the H an- 
tibonding orbitals which eventually leads to dissociation, 
and builds up on the molecule because the Mg surface is 
unwilling to accept back-donation of electrons from the 
H atoms. The Coulomb energy of this charge transfer is 
probably the main contribution for the energy barrier. 

By contrast, on the Ni doped surface there is almost 
no charge transfer from the substrate to the molecule at 

right). This is because 



the transition state (see Fig. 13 



while some Ni d charge fills up the H2 antibonding or- 
bitals, charge from the molecular bonding orbital is back- 
donated to the empty d states available on the surface. 
As a result, the energy barrier is reduced to almost zero. 



IV. EXPERIMENTAL 

A. Sample preparation 

The three samples prepared were MgH2, 2%Ni/MgH2 
and 2%Ti/MgH2. Three batches of 25 gram samples were 
prepared by ball milling the different compositions for 2 
hours under 4 bar of hydrogen. The MgH2 used for all is 
Goldschmidt 98% pure. The Ni used (99.9% pure) was 
from Alfa Aesar 0.8 - 0.03 /im diameter as was the Ti 
used (<20 fxm and 93% pure). 25 g of both 2%Ni/MgH 2 
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and 2%Ti/MgH2 were mixed in a tubular mixer before 
milling for one hour. Samples were then milled using the 
Fritsch planetary ball mill pulverisette 5. The milling 
pots have a special stainless steel jacket with an o-ring 
fitted on the top seal, this can allow a gas atmosphere 
through a feeding valve, to be used during the milling 
process of up to 5 bar. 25 grams of sample were milled 
using agate pots (around 300cc volume) and 15 balls of 
the same material. The milling process was 2 hours using 
350 rpm in a 15 minute mill 10 minute pause sequence. 



B. Sample testing 

The rig used for testing the sample has a 10 cc reactor 
pot containing 1 gram of sample. The main lines of the 
rig are a hydrogen line regulated to a 7 bar gauge, an 
Argon line and a vacuum line. The inert gas line and 
a vacuum line are used for purging the system. The re- 
actor is connected to an inlet flow controller, a pressure 
transducer to read the internal pressure and a mass flow 
meter outlet. A thermocouple in close contact with the 
powder load reads the sample temperature. A heating 
jacket cartridge is attached to the reactor allowing the 
system to operate in isothermal conditions or be tem- 
perature programmed from a control box which uses the 
sample thermocouple as a reference value. An interface 
card records inlet flow, outlet flow, temperature and pres- 
sure every second. 

Volumetric hydrogenation and dehydrogenation cycles 
were possible to monitor using this arrangement. Hydro- 
genations were performed at 300 C using 25 cc min -1 
of hydrogen from a regulator set at 7 bar gauge. During 
the hydrogenation the pressure increases to a point where 
the sample starts absorbing and forms a plateau pressure, 
once the sample is fully hydrogenated the reactor keeps 
building pressure until seven bar are reached. Dehydro- 
genations were recorded by flowing 25cc/min of hydrogen 
through the reactor using an inlet flow controller. The 
system is then open to vent for a chosen temperature or 
heating slope and any hydrogen evolving from the sam- 
ple is recorded as an increase in the 25cc min flow by an 
outlet flowmeter. 



C. Experimental Results 

The hydrogenation plots of the 2% Ni/Mg sample at 
different temperatures in the range between 290 and 
320 C are shown in Fig. [l4j The hydrogenation of the 
2%Ti/Mg sample in the temperature range between 290 

Both graphs of hydro- 



and 310 C is shown in Fig. 15 



genation show P mbar gauge vs time in seconds. The 
hydrogenation of pure Mg gave results close to the Ti 
doped samples and is therefore not shown. 

The catalytic activation of Mg by Ni during hydro- 
genation is clear. The plateau pressures for each temper- 
ature hydrogenation are lower for the 2%Ni/Mg sample 



than those of Mg and 2%Ti/Mg. The fact that at 290 C 
the hydrogenation curve for 2%Ni/Mg is still lower than 
the hydrogenation in the same conditions at 300 C con- 
firms this. At 290 C both Mg and 2%Ti/Mg showed a 
higher hydrogenation pressure than at 300 C suggesting 
Ti catalysisof Mg hydrogenation is not evident. 



V. CONCLUSIONS 

We have presented a DFT study of hydrogen dissocia- 
tion and diffusion over Ni-doped and Ti-doped Mg(0001) 
surfaces, and compared these with dissociation and diffu- 
sion on pure Mg(0001). Our results show that the energy 
barrier for hydrogen dissociation is high on the pure Mg 
surface (0.87 eV), and it is small (0.06 eV) or even null 
when Ni or Ti are used as dopants respectively. We also 
found that although on the pure Mg(0001) surface the 
binding energy of the H fragments is nearly zero, on the 
Ni/Ti -doped surfaces this binding energy is significant, 
being 0.66 and 1.34 eV respectively, with diffusion energy 
barriers of ~ 0.18, ~ 0.27 and ~ 0.8eV on the pure Mg, 
Ni-doped and Ti-doped surfaces respectively. 

Interestingly, the activation barriers for H2 dissociation 
over the Ni/Ti-doped Mg surface are similar to the values 
found on the corresponding pure Ni/Ti surfaces [551 1551 
|5J. 

More insight in the behaviour of these systems can be 
gained by inspecting the partial density of states and by 
looking at the electronic charge density distributions. In 
particular, the higher reactivity of Ti with respect to Ni 
can be understood in terms of a lower position of the ri- 
band centre, which correlates with both the height of the 
energy barriers for the dissociation of the H2 molecule 
and with the binding energy of the H fragments when 
adsorbed on the surface. 

The charge density distributions on the different sys- 
tems also shows some interesting behaviour. In partic- 
ular, we argued that the presence of a barrier on the 
pure Mg(0001) surface may be understood in terms of the 
build up of extra charge on the H2 molecule as it moves 
closer to the surface. This happens because the closed 
shell Mg surface is unwilling to accept back-donation of 
charge from the H 2 molecule. One consequence of this is 
that the molecule needs to arrive very close to the sur- 
face before starting to dissociate. By contrast, Ni and 
even more so Ti have many available empty d-states, 
and this avoids significant charge transfer from the sub- 
strate to the molecule. In this case, the dissociation of 
the molecule begins much further away from the surface. 

The low dissociation barrier, coupled with the low dif- 
fusion barrier, make Ni a very useful promoter for the 
hydrogenation of Mg. By contrast, the high dissociation 
barrier on the pure Mg surface, and the high diffusion 
barrier on the Ti-doped surface, are responsible for the 
slow kinetics of hydrogenation on both systems. 

Our experimental findings show faster hydrogenation 
for the 2%Ni/Mg sample with respect to the reference Mg 
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or the 2%Ti/Mg, in good agreement with our theoretical 
results of a lower activation energy for the dissociation- 
diffusion process in the 2%Ni/Mg system. The behaviour 
of Mg and 2%Ti/Mg upon hydrogenation is found to be 
very similar, again agreeing very well with the theoretical 
findings of large and similar activation energies: a disso- 
ciation energy barrier of 0.87 in the pure Mg system, and 
a diffusion energy barrier of 0.8 eV in the 2%Ti/Mg sys- 
tem, making the dissociation-diffusion process similarly 
difficult in both cases. 

We deliberately chose to study Ni and Ti as dopants 
because they are at the two ends of the first row of tran- 



sition metals, and so their behaviour may be expected 
to be representative of a range of properties. In fact, we 
are now extending our investigations to other transition 
metals, and we plan to report on these new results in the 
near future. 
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FIG. 1: (Colour) Possible adsorption sites (top, bridge, 
hollow-hep and hollow-fcc) for hydrogen (dark red) on the 
Mg(0001) surface (light blue). 
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FIG. 2: Minimum Energy Path for H2 dissociation and dif- 
fusion on a pure Mg(0001), Ni-doped Mg(0001) and Ti-doped 
Mg(0001) surface. 



FIG. 3: (Colour) H2 (dark red) dissociation on the pure 
Mg (light blue) surface as viewed from side (top figures) and 
top (bottom figures). Figures show positions at IS (left-hand 
panel), TS (central panel) and FS (right-hand panel). 




FIG. 4: (Colour) H (dark red) diffusion on the pure Mg (light 
blue) surface as viewed from top. Figures show positions at 
FS (top-left), TS2 (top-centre), LS (top-right), TS3 (bottom- 
left) and FS2 (bottom-right). 
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FIG. 5: (Colour) Possible final state adsorption sites for H2 
(dark red) dissociation over the metal-doped (Ti/Ni) (dark 
green) Mg surface (light blue). In the case of the Ni doped 
surface the bottom-left site was not a stable configuration. 




FIG. 6: (Colour) Same as Fig. [3] but for H2 dissociating over 
the Ti-doped Mg surface at IS and FS (there is no TS in this 
case). The Mg, Ti and H atoms are represented respectively 
by light blue, dark green and dark red colours. 



FIG. 7: (Colour) Same as Fig. [i] but for H diffusion over 
the Ti-doped Mg surface. Figures show positions at FS (left), 
TS2 (centre) and FS2 (right). The Mg, Ti and H atoms are 
represented respectively by light blue, dark green and dark 
red colours. 




FIG. 8: (Colour) Same as Fig. [3] but for H2 dissociating 
over the Ni-doped Mg surface. The Mg, Ni and H atoms are 
represented respectively by light blue, dark green and dark 
red colours. 



FIG. 9: (Colour) Same as Fig. [I] but for H diffusion over 
the Ni-doped Mg surface. Figures show positions at FS (left), 
TS2 (centre) and FS2 (right). The Mg, Ni and H atoms are 
represented respectively by light blue, dark green and dark 
red colours. 




FIG. 10: Projected density of states for H2 dissociating over 
a pure Mg surface as a function of the energy relative to the 
Fermi level, respectively for the initial state (IS, top-left cor- 
ner), transition state (TS; top-right corner), transition state 
plus one further step along the MEP (TS+1; bottom-left cor- 
ner) and final state (FS; bottom-right corner). 
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FIG. 11: As in Fig. [lO]but for the Ni-doped Mg surface. The 
dashed vertical line shows the position of the d-band centre. 
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FIG. 13: (Colour) Charge distribution during H2 (dark red) 
dissociation at the TS of the MEP respectively on the pure 
Mg (left) and Ni-doped Mg (right) surfaces (see text for de- 
tails). White shows positive charge and black negative charge. 
Isolines are also shown in white. 





1,1,1 




k 


6000 




— 290 C 
■ ■ ■ 300 C 
■- 310 C 

- 320 C 


14 

// 


Is 4000 

S 

a. 






ft 


2000 


" f> "* 












A 









\) 500 1000 1500 

Time (sec) 



FIG. 14: (Colour) Different temperature hydrogenation plots 
for 1 gram of 2%Ni/Mg using 25 cc/min of H2. 
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FIG. 15: (Colour) Different temperature hydrogenation plots 
for 1 gram of 2%Ti/Mg using 25 cc/min of H2. 



